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Abstract 

The Einstein equations for one of the hypersurface-homogeneous rotating dust mod- 
els are investigated. It is a Bianchi type V model in which one of the Killing fields is 
spanned on velocity and rotation (case 1.2.2.2 in the classification scheme of the earlier 
papers). A first integral of the field equations is found, and with a special value of 
this integral coordinate transformations are used to eliminate two components of the 
metric. The k = — 1 Friedmann model is shown to be contained among the solutions in 
the limit of zero rotation. The field equations for the simplified metric are reduced to 3 
second-order ordinary differential equations that determine 3 metric components plus 
a first integral that algebraically determines the fourth component. First derivatives 
of the metric components are subject to a constraint (a second-degree polynomial with 
coefficients depending on the functions). It is shown that the set does not follow from 
a Lagrangian of the Hilbert type. The group of Lie point-symmetries of the set is 
found, it is two-dimensional noncommutative. Finally a method of searching for first 
integrals (for sets of differential equations) that are polynomials of degree 1 or 2 in the 
first derivatives is applied. No such first integrals exist. The method is used to find a 
constraint (of degree 1 in first derivatives) that could be imposed on the metric, but it 
leads to a vacuum solution, and so is of no interest for cosmology. 



1 Statement of the problem and summary of the paper. 

This paper is a continuation of a series of papers on rotating dust models in relativity 1-3 . 
The initial motivation for this research was the desire to find a rotating generalization of 
the Friedmann models. In spite of much effort spent on investigating solutions of Einstein 
equations with a rotating matter source, no such generalization has been found so far; see 
literature surveys in Refs. 3 and 4. Refs. 1, 2 and 3 provided a complete classification 
scheme for hypersurface-homogeneous rotating perfect fluid models with zero acceleration. 
Unlike in previous approaches, nothing was assumed about the position of the symmetry 
orbits in spacetime; the classification includes also timelike and null orbits, and so it is the 
farthest-reaching application of the Bianchi classification to rotating and nonaccelerating 
perfect fluid models in relativity. The models split into 3 general classes: I, in which two 
of the Killing fields are everywhere spanned on the vector fields of velocity u a and rotation 
ii; a (Ref. 1); II, in which only one Killing field is spanned on u a and u> a (Ref. 2); and 
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Ill, in which all Killing fields are linearly independent of u a and w a (Re£. 3). The many 
particular cases arise because of several possible alignments or misalignments among the 3 
Killing fields and u a and w a . 

By the Bianchi type of the symmetry algebra and by the relation of the velocity field to 
the symmetry orbits it can be recognized in which cases generalizations of the Friedmann 
models can be expected. Two such candidate cases were found in class II, and five more 
in class III. Those of class III were prohibitively complicated, but one of the cases of class 
II allowed for some progress, and this one is presented in the present paper. It is the 
Bianchi type V subcase of the case 1.2.2.2, given by eq. (5.19) in Ref. 2. The other 
candidate case found in class II, eq. (5.10) in Ref. 2, can reproduce only the de Sitter or 
the Einstein model in the limit of zero rotation, this is seen from the time-dependence of 
the metric. Hence, it is not interesting for cosmology and therefore disregarded here. In 
sec. 2, the metric is simplified by a coordinate transformation, and a first integral of the 
Einstein equations is found. With zero value of this integral, coordinate transformations 
can be used to eliminate two components of the metric tensor, and the number of nontrivial 
Einstein equations is reduced to 7. Although there are only 4 functions + matter density to 
be determined by these 7 equations, the set later turns out to be self-consistent. In sec. 3, it 
is shown that the k = — 1 Friedmann models are contained among the solutions that result 
in the limit of zero rotation. In sec. 4, the Einstein equations are reduced to a set S of 3 
second-order equations to determine 3 metric components + a quadrature Q to determine 
the fourth component (533). Of the Einstein equations derived in sec. 2, one is fulfilled 
identically in consequence of the set {S{J Q}, one turns out to be a constraint imposed on 
the initial data, and the one that determines the matter-density turns out to provide a first 
integral. The constraint and the first integral are second-degree polynomials in the first 
derivatives of the unknown functions whose coefficients depend on the unknown functions. 
The first integral determines 533 algebraically in terms of the other components, and so it 
is a replacement for the quadrature Q. It is also shown that the set S cannot be obtained 
as the Euler-Lagrange equations from a variational principle of the Hilbert type. Finally, it 
is shown in sec. 4 how the set {S[jQ} reproduces the Friedmann equations in the limit of 
zero rotation and zero shear. In sec. 5, Lie point-symmetries of the set are found: there is 
a two-dimensional symmetry group that allows one to reduce one second-order equation to 
a first-order equation plus a quadrature. However, this reduction provides no real progress 
toward solving the set S; the first-order equation is still a member of a complicated set. 
In sec. 6, a method of systematic search for polynomial first-order first integrals of a set 
of ordinary differential equations is applied to the set S of sec. 4. It is shown that no 
first integrals that are polynomials of degree 1 or 2 in the first derivatives exist. The same 
method is used to reveal the existence of a possible constraint on initial data, which is 
of degree 1 in first derivatives, that is preserved by the set S. However, the constraint 
necessarily implies zero matter-density, and so it is not interesting for cosmology. 

Calculations that are of secondary importance for the main text, but are difficult to 
reproduce, are described in the appendices. 
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2 The Einstein equations, their first integral and implica- 
tions of the zero value of this integral. 

The subject of the present paper are the Einstein equations for the Bianchi type V subcase 
of case 1.2.2.2 of Ref. 2. For reference, the initial formulae are recalled in their original 
notation. 

The Bianchi type V symmetry results when c = in eqs. (5.19) of Ref. 2 and when, in 
addition, j = —a in eqs. (5.16). Hence, the metric is: 

ds 2 = di 2 + 2ydtdx + y 2 h u dx 2 + 2h 12 dxdy + 2y 2 h l3 dxdz 

+ (h 22 /y 2 )dy 2 + 2h 23 dydz + y 2 h 33 dz 2 , (2.1) 

where the coordinates are {x a } = {x°, x 1 , x 2 , x 3 } = {t,x,y,z}, and hij,i,j = 1,2,3 are 
unknown functions of the variable 

v = e*2/ C2/a , (2.2) 

a and C 2 being arbitrary constants. The velocity field u a , the rotation field w a and the 
Killing fields k^ a , i = 1, 2, 3 are given by: 

u a = 6 a , w a = (p/y)S a , k (1 f = 6 a 1 , k (3) a = 5 a 3 , 

k {2) a = C 2 5 a + a(x5 a 1 - yb a 2 + zS a 3 ), (2.3) 

where p is the matter-density of dust. The rotation tensor uj a p has only one algebraically 
independent nonzero component: 

^ = \ ( 2 - 4 ) 

and therefore the coordinates used here are ill-suited for considering the limit to — > 0. From 
the first equation in (2.3) it can be seen that the coordinates are comoving. 

As shown in Ref. 1, it follows from the equations of motion and from the equation of 
conservation of the number of particles that: 

g := det( 5a/3 ) = -(y/p) 2 , (2.5) 

where p is the mass-density. 

This is the form in which the metric resulted from the Killing equations in Ref. 2. It 
is advantageous to transform the coordinates as follows: 

t = t'-(C 2 /a)lny', x = x' - C 2 /(ay'), (y,z) = (y',z'). (2.6) 

The result is equivalent to substituting C 2 = and a = 1 in eqs. (2.1) - (2.4), i.e. the 
forms of the metric (2.1), of the vector fields u a , w a , and k^ a in (2.3) and of the 

rotation tensor u> a p in (2.4) do not change (although the new h'^ in (2.1) will be linear 
combinations of the old hij), while the new k^ a basis vector will be: 

k {2 f = x5 a l -y5 a 2 + z5 a 3 , (2.7) 
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and the argument of will now be v = e*' , i.e. the hij are from now on unknown functions 
of the time-coordinate t. 

The isometry corresponding to (2.7) is: 

t' = t, (x',z') = e T (x,z), y' = e- T y, (2.8) 

where r is the group parameter. 

It is convenient to parametrize the metric as follows: 



ds 2 = {dt + ydx) 2 -{yK^dx) 2 -{K/y) 2 {dy + y 2 hdx) 2 -K 33 2 [ygdx+{f /y)dy + ydz] 2 , (2.9) 

where Ku, K, K 33 , h, /, and g are unknown functions of t. The components of the Einstein 
tensor referred to below are tetrad components Gjj = e a ie 13 ,jG a p, i.e. projections of the 
coordinate components G a p onto the orthonormal tetrad e 1 := e 1 a dx a implied by (2.9): 

e° = dt + ydx, e 1 = yKudx, e 2 = (K/y)(dy + y 2 hdx), 



e 3 = K 33 [ygdx + (f/y)dy + ydz], (2.10) 

where e a j is the inverse matrix to e 1 a , i.e. e a je I p = 5 a p, e a je I a = S 1 j. In the parametriza- 
tion (2.9), the determinant of the metric is: 

g=-(yK ll KK 33 ) 2 . (2.11) 

The tetrad components of the Einstein tensor corresponding to the metric (2.9) are 
given in the Appendix A. As seen from there, two combinations of those equations are of 
first order, they are i^nGo3 + G\ 3 = 0, i.e.: 

(^33/^ii)[(^n 2 - l)K' 2 f, t +h(hf, t -g, t )} = (2.12) 
and KuGo2 + G12 = 0, i.e.: 



{K ll K)- 1 [~K 2 hh, t +h - K n K n ,t +(K U 2 - l)(2K, t /K - K 33 , t /K 33 )\ = 0. (2.13) 

As shown in Appendix B, the case h = does not lead to interesting developments, so we 
shall proceed further under the assumption: 



h^O. (2.14) 

Then, eq. (2.12) implies: 

g, t =[h + (K n 2 -l)/(hK 2 )}f, t (2.15) 

With this, the equations G03 = G13 = G23 = turn out to be equivalent, and they can be 
written as follows: 

1 / K u 2 - 1 i^33 3 /,A K 33 3 f, t , 01 m 

" 2 { —fT- • ~K^K ) " + 1^K ~ °" (2 ' 16) 
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This invites the introduction of the new variable u(t) by u, t = h/ (Kn 2 — 1), and then (2.16) 
becomes: 

K U K ) K llK ~ °' (A17j 

which has the first integral K 33 3 /,« /(KnK) = Ce 2u , C = const, i.e.: 

f,t = Ce 2u hKnK/[K 33 3 (Kn 2 - 1)]. (2.18) 

From here on, we shall follow only the special case C = 0, which is a solution of the Einstein 
equations, but not a general one: it is a subcase chosen ad hoc for further progress with 
integration. Then, from (2.18) and (2.15) / = const, g = const, and from (2.10) the 
coordinate transformation z' = z + / / y + gx leads to 

f = g = (2.19) 

without changing any of the other formulae for g a p, u a , w a , u) a p or k^ a . 

The Einstein equations G03 = G13 = G23 = are now fulfilled identically. We are 
left with 7 equations of the set (A.l) - (A. 10) in Appendix A that should determine the 
4 functions Kn, K, K33 and h, and the matter density p in addition. It will turn out in 
sec. 4 that the 7 equations are dependent just in the way needed to make the problem 
self-consistent and determinate. 



3 The Friedmann limit of the metric. 

As already stated, the coordinates used in sec. 2 are ill-suited for considering the limit u — > 
0. It will be shown in the present section that this limit can be calculated after a coordinate 
transformation and a reparametrization of the metric. This is just a demonstration of 
existence, and it is not claimed that the limit oj — > thus obtained is unique (i.e. another 
nonrotating limit might be obtained starting from a different coordinate transformation). 
However, we will be satisfied to show that a limit exists in which the k < Friedmann 
model is contained. 

Since U12 = —OJ21 = \ are the only nonzero components of the rotation tensor, a natural 
coordinate transformation to consider is: 

V = u y'. (3.1) 
where u>q is a constant. After the transformation: 

w i2 = ^0 = ( 3 - 2 ) 

(all other u a p = 0), and the limit of zero rotation is ujq — ► 0. However, before this limit 
is taken, the metric functions in (2.9) must be reparametrized or else the limit will be 
singular. The following reparametrizations will do the job: 

K u = Ku/uq, K 33 = K 33 /uj , f = fu . (3.3) 
The transformation (3.1) and the reparatmetrization (3.3) result in the following metric: 
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ds 2 = (dt + u y'dx) 2 - (y'Kndx) 2 - K 2 (dy'/y' + uJ Q y'hdxf 
-K 2 3 [y'gdx + (f/y')dy' + y'dz} 2 
whose limit — > (with primes and tildes omitted) is: 



(3.4) 



ds 2 = dt 2 - (yK u dx) 2 - {K/yfdy 2 - K 33 2 [ygdx + (f/y)dy + ydz} 2 , 



(3.5) 



This is more than sufficiently general to accomodate the k = — 1 Friedmann model that 
results when g = f = and K\\ = K = K33 := R(t), where R(t) is the Friedmann scale 
factor. The resulting coordinates are none of the standard ones, but are related by y = e u 
to one of the sets used in the literature (see eq. (1.3.15) in Ref. 5). 

The fact that (3.5), the limit luq — > of (2.9), is still more general than the Friedmann 
metric means that (3.5) has nonzero shear, i.e. shear survives the transition uj — > 0. 

However, one possible problem still lies ahead. It was proven above that the k = — 1 
Friedmann model is contained among the solutions of the set (A.l) - (A. 10). What is still 
needed is an explicit solution with the property that it has nonzero rotation in general, 
but reproduces the k = — 1 Friedmann model in the limit to — > 0. Experience with the 
Einstein equations in other cases shows that sometimes, while integrating the equations, 
one encounters mutually exclusive alternatives A and B such that it is no longer possible 
to recover B as a limit of A after the integration is completed. A well-known example are 
the two subfamilies (/?' = and f3' 7^ 0) of the Szekeres-Szafron 6 ^ 7 cosmological models; 
see Ref. 5 for more on this point. (Only recently, a reformulation of the two classes 
was invented that allows to recover the 0' = family from the other one, see Ref. 8). 
Hence, it may still happen that among the explicit solutions, the rotating dust model and 
the Friedmann k = — 1 model will turn out to be mutually exclusive subfamilies. This 
uncertainty will persist until an explicit solution is found. 

It will be shown at the end of sec. 4 that the explicitly written out Einstein equations 
do allow a continuous limiting transition uj — > 0, a — > 0, and in the limit they reproduce 
exactly the Friedmann equations. 

4 The independent Einstein equations. 

We shall now proceed with the subcase (2.19). Eq. (2.12) is then fulfilled identically. Eq. 
(2.13) does not change, and it can be more conveniently rewritten if K\\ is parametrized 
as follows: 



When this is substituted into the remaining equations (A.l) - (A. 10), the function K33 
disappears from the set completely, i.e. we are left with 6 equations to determine h, K, F 
and the matter-density plus the quadrature implied by (4.2) that allows one to calculate 
If 33 once h(t), K(t) and F(t) are known. 



K u = cosh(F). 



(4.1) 



Then, from (2.13): 
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Since (2.13) is now satisfied, the equations Gq2 = and G±2 = are equivalent, and 
they can be written as: 

h, tt = ^K 2 hh, t 2 /smh 2 (F) - 5K, t h, t /K + (2cosh 2 (F) - l)F, t h, t /sinh(F) cosh(F) 



+hh, t /sinh 2 (F) + K, t /K 3 + F, t /K 2 cosh(F)sinh(F) - -h/(Ksmh(F)) 2 (4.3) 

This is used to eliminate h,tt from the other Einstein equations. The equation Goi = can 
then be solved for F,tt (the solution is given in Appendix C) and this is used to eliminate 
F,tt from the diagonal components of the Einstein tensor (all the non-diagonal Einstein 
equations have been used up at this point). After such a substitution, the following identity 
is fulfilled: 

Gn + G 33 - 2G 22 = 0, (4.4) 

i.e. one of the three equations Gn = G22 = G33 = A can be discarded because it is a 
consequence of the remaining two. We choose to discard G33 = A. 
Then, K,u can be calculated from G22 — Gn = 0. The result is: 

K,tt=^K 3 smh- 2 (F)h, t 2 - ^K 3 hcosh(F)smh- 3 (F)F, t h, t 

-cosh- 1 (F)smh'\F)K, t F, t +2 cosh(F)sinh- 1 (F)K, t F, t -Kcosh 2 (F)smh- 2 (F)F, t 2 
~Kh, t +^K 3 h 2 smh- 4 (F)h, t +^K 3 h 2 smh- 2 {F)h, t 

~hK, t -hsmh.- 2 (F)K, t -Khcosh- 1 (F)smh.- 1 (F)F, t +^Khcosh 3 (F)smh- 3 (F)F, t 

--if _1 cosh 2 (F)sinh- 2 (F) - -Kh 2 sm\T i {F) - ^/i 2 sinh" 2 (F) (4.5) 

This is used to eliminate K,u from the right-hand side of the equation determinig F,tt (see 
Appendix C), and the result is: 

F, tt = -^K 2 cos\r 1 (F)s\nh- 1 {F)h, t 2 - ^Khcosk' 1 {F)s\n\r 1 {F)K, t h, t 
+2K- 2 cosh~ 1 (F) sinh(F)K, t 2 - K~ x K, t F, t - cosh(F)sinh~ 1 (F)F, t 2 

+^ 2 /i 2 cosh~ 1 ( J F)sinh" 1 (F)/i, t +^ 2 /i 2 cosh- 1 ( J F)sinh- 3 (F)/i, t 

1 3 

+-cosh _1 (F)sinh _1 (F)ft,,j --cosh -1 (.F) smh(F)h, t 
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~K- 1 hcoBh.- 1 (F)gmir 1 (F)K, t —K^hcosh^iF) smh(F)K, t 



+hsmh- 2 (F)F, t +^hF, t -^- 2 cosh _1 (F) sinh(F) - ^K~ 2 cosh- 1 (F)sinh- 1 (F) 

-^ 2 cosh- 1 ( J F)sinh _3 (F) - i/i 2 cosh~ 1 (F)sinh~ 1 (F) (4.6) 

With (4.2), (4.3), (4.5) and (4.6) all substituted into (A. 5), the equation Gu = A reduces 
to the following form: 

Gu = ^K 2 cosh- 2 (F)/i,i 2 + ^Khcosh~ 2 (F)K, t h, t +^K 2 hcosh- 1 (F)smh- 1 (F)F lt h, t 
-2K- 2 cosh- 2 (F)smh 2 (F)K, t 2 - 2J^ 1 cosh- 1 (F) smh{F)F, t K, t +F, t 2 
+^K 2 h 2 cosh~ 2 (F)h, t -3K 2 h 2 smh- 2 (F)h, t -^cosh" 2 (F)h, t +^h, t 



+^K- 1 hcosh- 2 (F)K, t +3K~ l hK, t ~hcoBh- 1 {F)snh.- 1 (F)F, t -3fccosh _1 (.F) gmh(F)F, t 

+^K- 2 cosh- 2 (F) + ^K~ 2 + ^h 2 cosh- 2 (F) + /i 2 sinh~ 2 (F) = A. (4.7) 

Now it may be verified that Gu = const is preserved by eqs. (4.3), (4.5) and (4.6). This 
is done as follows. The derivative ^Gn is calculated, and h, tt , K, tt and F tt that reappear 
are eliminated using (4.3), (4.5) and (4.6). Then, K, t 2 is found from (4.7) and used to 
eliminate K,t 3 and K, 2 from ^Gu. The result is the identity ^Gn = 0. This means 
that, in virtue of the other field equations, if Gu = A holds at any given time, then it will 
remain constant at all other times. Hence, Gu = A is a limitation imposed by the Einstein 
equations on the initial data for eqs. {(4.3), (4.5), (4.6)}, and it defines the cosmological 
constant in terms of the other constants that will appear after (4.3), (4.5) and (4.6) are 
solved. If A = 0, then Gu = reduces the number of arbitrary constants by 1. 

Hence, with (4.4), we are left with only four equations: (4.2), (4.7) and any two equa- 
tions from the set S = {(4.3), (4.5), (4.6)}, to determine the four functions K33, h, K and 
F. The third equation in S is implied by the remaining two together with (4.7). The only 
field equation that has not yet been used up is: 

Goo = (87rG/c 4 )p - A. (4.8) 

This may be expected to simply define the matter-density in terms of the metric functions. 
However, in the formulation used in this paper, matter-density enters the equations in two 
ways: as a source term in Goo above, and also through (2.5). From (2.5) and (2.11) it 
follows that p must be related to the other functions by: 

p = (K 11 KK 33 )- 1 . (4.9) 
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Together with (4.8) and (4.1) this implies that the following must hold: 

[(Goo + A) cosh(F) KK 33 ], t = 0. (4.10) 

Indeed, this is an identity. This is verified as follows. First, (4.1), (4.2), (4.3), (4.5) and 
(4.6) are substituted into (A.l) (with / = g = 0) to eliminate all second derivatives. Then, 
(4.10) is calculated, and (4.2), (4.3), (4.5) and (4.6) are used to eliminate K 33 ,t and all 
second derivatives again. Finally, (4.7) is used to eliminate K,t S and K, 2 from the left-hand 
side of (4.10). In the end, the identity (4.10) results. Hence, (2.5) and (4.8) are consistent 
in virtue of the other field equations, and moreover (Goo + A) cosh(F)KK 3 3 = G = const 
(with second derivatives of h, K and F eliminated by (4.3), (4.5) and (4.6) and with K, 2 
eliminated by (4.7)) is the following first integral of the Einstein equations: 

K 33 [-3Khsmh(F)F, t --Kh 2 cosh~ l (F) + ^Kcosh- 1 (F)smh 2 {F)h, t 

+3/icosh" 1 (F)sinh 2 ( J F)i^,i -^- 1 sinh 2 ( J F)cosh~ 1 (F) - -K 3 h 2 cosh' 1 (F)h, t ] = C. 

(4.11) 

Note that, from (4.8) and (4.9), G = SttG/c 4 ^ 0, and so (4.11) determines K 33 alge- 
braically. Hence, (4.11) can replace (4.2) as the definition of ^33. Thereby, the problem of 
this paper was reduced to the following procedure: 

1. Find the most general solution of the set {(4.3), (4.5), (4.6)}. It will contain 6 
arbitrary constants {Gi, . . . , Cq}. 

2. Impose (4.7) on the {h, K, F} found in the previous step. This will be just a 
definition of A in terms of {Gi, . . . , Cq} or, when A = 0, an additional constraint imposed 
on {Gi, . . . ,G 6 }. 

3. Calculate K 33 from (4.11), with G = 8ttG/c 4 . 

4. Calculate the matter-density from (4.9). 

As shown in Ref. 9, an efficient method to find first integrals of a set of equations exists 
if the set can be obtained from a Lagrangian. Unfortunately, the problem of determinig 
whether a given set of equations is derivable from a lagrangian is rather complicated and 
unsolved in general 10 . It is known that the Einstein equations for class B Bianchi metrics 
may not admit a lagrangian, even though the general Einstein equations do (see Ref. 11 
for an explanation). It is shown in Appendix D that eqs. {4.3), (4.5), (4.6)} do not follow 
from the most natural lagrangian conceivable in this CclSCI Si second-degree polynomial in 
the first derivatives of h, K and F, with coefficients being functions of h, K and F. 

For further reference, let us consider the limit of zero rotation in (4.2) - (4.3) and (4.5) 
- (4.7). After the reparametrization (3.3) we have: 

cosh(F) = Ku/uo, sinh(F) = \f K n %o 2 - 1, 

F, t = iWV^n 2 -^) 2 , (4.12) 
and then (4.2) in the limit loq — > becomes: 

K 33>t = K 33 (2K, t IK - K 11>t /K u ), (4.13) 
which is an identity in the Friedmann limit Kn = K = K 33 = R(t). 
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Eq. (4.3) could in fact be discarded in the limit uq — > 0. This is because eq. (4.3) was 
derived from (A. 3), and those terms in (A. 3) that lead to (4.3) are all multiplied by ojq 2 
after the reparametrization (3.3). The off-diagonal component of (3.4) that is proportional 
to h will vanish with any h when ujq — > 0. Nevertheless, (4.3) gives a result consistent with 
the other equations in this limit. The limiting form of it is: 

h, tt = -5h, t K, t /K + 2K lu h, t /K~ u + K, t /K 3 . (4.14) 
The limit lo -> of (4.5) is: 

q Q Q 

K, tt = -KK llt t /Kn 2 ~ -Kh,t--hK, t +-KhK lltt /K 11 - 1/(4K) + 2K, t K n , t /K n . 

(4.15) 

The same limit of (4.6) is: 

Kn,tt/K u = 2K, t 2 /K 2 - K, t K n , t /(KK n ) + ^hK lltt /K u - h, t ~hK, t /K - l/(4K 2 ). 

(4.16) 

In the Friedmann limit Kn = K = R(t), eqs. (4.15) and (4.16) become identical: 

R, tt /R = R, t 2 /R 2 - ^h, t -1/(4R 2 ). (4.17) 
Finally, the limit loq — > of (4.7) is: 

- 2K, t 2 /K 2 - 2K, t Kn,t/(KK n ) + Kn/ / K u 2 + -h, t +3hK, t /K - 3hK n ,t/K n 

+3/(2K 2 ) = A. (4.18) 

The Friedmann limit of this is: 

-3R, t 2 /R 2 + h, t +3/(2R 2 )=A. (4.19) 
Finding h, t from (4.19) and substituting it in (4.17) we obtain: 

R,u /R = -R,t 2 /(2R 2 ) + l/(2R 2 ) - A/2, (4.20) 

which is exactly one of the Friedmann equations. Incidentally, the h,t found from (4.19), if 
substituted in (4.14), leads to (4.20) again. Hence, in the Friedmann limit, (4.14) follows 
from (4.19) and (4.17), and need not be discarded. 

Note that also (4.11) has a meaningful Friedmann limit. In order to make this limit 
finite, it must be assumed that: 

C = C/uo 2 , (4.21) 

and then the limit ujq — ► of (4.11) is: 



K 33 [AKK U + 2K, t 2 K u /K + 2K, t K 11;t - KK u ,t /K n - 3K U /K] = C. (4.22) 
In the Friedmann limit this becomes: 
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R(AR 2 + 3R, t 2 -3) = C. 



(4.23) 



Recalling the Friedmann formula for the mass-density, with k = — 1: 

3R, t 2 /R 2 - 3/R 2 + A = (8irG/c 2 )p, (4.24) 

we recognize in (4.23) the familiar mass-conservation formula of the Friedmann model, 
pR? = c 2 C/(8itG) = const. 



5 The Lie point-symmetries of the equations (4.3), (4.5) and 
(4.6). 

Point symmetries of (sets of) differential equations are transformations in the space of 
the independent + dependent variables that leave the set of solutions of the equations 
unchanged. The point symmetries that form Lie groups (if they exist for a given set of 
equations) can help in transforming apparently intractable equations into solvable ones 
by adapting the variables suitably to the generators of the symmetries. The background 
philosophy and many of the methods are analogous to simplifying the Einstein equations 
by adapting the coordinates to the Killing vector fields (if such exist). It is assumed that 
the readers are familiar with this latter procedure. The basic definitions and theorems 
concerning point symmetries are presented in detail in Refs. 9 and 10. 
Eqs. (4.3), (4.5) and (4.6) are of the following form: 



dV 



• dz J dz k ■ dz- 7 



dt 2 _ " ]k dt dt 1 r 1 dt 1 ~ ' (5,1) 

where i = 0, 1, 2; (z°, z 1 , z 2 ) = (h, K, F) and W % jk, V l j and U l are functions of the z l , but 
not of t. (Incidentally, the independence of t of all these coefficients immediately implies 
one group of symmetries, t — > t' = t + s, where s is the group parameter. This group will 
emerge from the calculation below.) Let the following be a one-dimensional group of point 
transformations: 



t' = t'(t,{zi},r), z H = z H (t,{zJ},r), (5.2) 

where r is the group parameter and t = tq corresponds to the identity (so that t'(t, {z J }, tq) = 
t, etc.). The generators of this group (the field of vectors tangent to the orbits of the group 
(5.2)) are then: 



X 



,d_, j_d_ 



where: 



(5.3) 





d 


' a ' 




" d7 


z'i 



T=T0 

k. 



(5.4) 



(fc) 



The generator X is extended to arbitrary derivatives := z by the recursive formulae: 

d k zi d£ 



(o)i 

V =rr>, 



(k-l)j 

(k)j d V 
V = 



dt 



dt k dt' 



(5.5) 
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and by: 



(*) d , d 0)i d d , , 

a z a z 



The derivatives 4z in (5.5) are total derivatives, i.e. 



dt 

-/(t,{z},{z},...,{z}) = - + — + E z 

p=i a z 

and the order n to which the generator X has to be extended is equal to the highest order 
of derivatives in the set (5.1) (n = 2 in our case). A generator of a point symmetry obeys 
then: 

(n-l)i 

■ d r? df 

(n)i 

where 0* is the right-hand side of (5.1). (The right-hand side of (5.7) is the V as given 

by (5.5), but with ^£ replaced by f2 l from (5.1)). Eqs. (5.7) must be identities in all the 

(l)i (n-l)i 

derivatives z , . . . , z , and so they imply several separate equations to be obeyed by 
the £ and rf. 

The procedure in finding and exploiting point symmetries is thus the following: 

1. Find the general solution of (5.7) for X. Since the generators form a Lie algebra 
(see Ref. 9), the most general X will be spanned on a finite number of basis vector fields 

x ik) . 

2. Read off the basis Xn^ from that solution. 

3. Adapt the variables {t'(t, {z J }), z H (t, {z^})} to the basis fields Xr k \ so as to maximally 
simplify the equations. 

For our equations (5.1), eqs. (5.7) imply the following four relations: 

Ul +W j M C, 3 = 0, (5.8) 



V\kl = W^sT," + 2W i s{l r, s , k) -W s kl r,\ s +<5\^ s fe) £, s +V\ (l £, fc) +2^^), (5.9) 
where parentheses on indices denote symmetrization, 

= W\,— + -V\ ,ri s + -V\ri s , k --V s k n\, 
dtdz k s dt 2 ' 2 1 2 ' ' 



1 d£ 1 1 ■ d 2 f 

+^V\ f t + IPU +-8> k U°Z, 8 +-fk-Q$, ( 5 - 10 ) 



d ~w = yls % + u ^ s 7/8 - USr]l > s +2W % (5 - n) 



12 



The general solution of these equations (with W 1 m, V l k and U l read off from (4.3), (4.5) 
and (4.6)) is: 

where A and B are arbitrary constants. The proof that this is the most general solution is 
laborious but straightforward, it is given in Appendix E. Hence, our set of equations has a 
two-dimensional symmetry group whose generators are: 

and the corresponding finite symmetry transformations are: 

t' = t + n, (h',K',F') = (h,K,F); 

t' = e T H, h' = e~ T2 h, K' = e T2 K, F' = F, (5.14) 

where t\ and T2 are the group parameters. The first symmetry was self-evident, as already 
mentioned, and the second one can be verified by inspection of the equations (4.3), (4.5) 
and (4.6). 

Unfortunately, these symmetries do not lead to any discernible simplification of the 
set S = {(4.3), (4.5), (4.6)}. In variables adapted to the generator the independent 

variable is K, and t(K) is one of the functions. The set (5.1) thus transformed is of first 
order in (f)(K) := dt/dK, but the first-order equation is still a member of a complicated set 
and none of the equations separates out. Moreover, after the transformed set is algebraically 
solved for t,KK, h,KK and F,kk, the right-hand sides become polynomials of third degree 
in t, k , h, K and F, K . 

The variables adapted to the generator X^ 2 ) are (t',h',K'), where: 



t', K = e K \ h = e- K 'ti. (5.15) 



In these variables, the set (5.1) becomes of first order in ip(t') = K', t '- However, after it is 
solved for h', t 't'i K' \ff and F, t > t >, the right-hand sides of h', t 't' and F, t > t > contain rational 
functions of the form W/(l + t'K', t i ), where W is a monomial of second degree in some 
of the h',f, K', t > and F, t i. Neither equation separates out. It is not possible to adapt 
the variables to both the generators simultaneously because the group is nonabelian. This 
author was not able to make any use of the new variables. 

6 First integrals that are polynomials in (h, t , K, t , F, t ). 

Suppose that the set S = {(4.3), (4.5), (4.6), (4.7)} has a first integral of the form: 

/ := QijZ^i + Li? + E = C = const, (6.1) 

where C is an arbitrary constant, Qij = Qji, L; L and E are unknown functions of (h, K, F), 
i,j = 1,2,3, z 1 = h,z 2 = K,z 3 = F. Then ^ = in virtue of S, i.e. using (5.1) to 
eliminate z l : 
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(2Q ij z> + Li)(W l kl z k z l + V\z k + U l ) + Q ijtk z i z j z k + L^z l z? + E» z { = 0. (6.2) 



In showing that (6.2) is zero, (4.7) must be used. Eq. (4.7) may be safely used to eliminate 
F, t 3 and F, t 2 , but not the remaining F, t . This is because F, t found from (4.7) would be of 
the form: 



where P and A are polynomials of degree 2 in h,t and K,t- If A were a square of a 
first-degree polynomial, then (6.3) could be used to eliminate F, t altogether from (6.2). 
However, A being a square implies an additional equation obeyed by h, t and K, t (the 
discriminant of A must be zero). Hence, if h,t and K,t are to be treated as independent, 
then F,t is linearly independent of h,t and K,t- Then the coefficients of F,t in (6.2) must 
sum up to zero anyway, and eliminating F, t is of no use. 

Knowing this, it can be verified that first integrals of the form (6.1) do not exist. 
The calculations are conceptually straightforward, but lead through horrible intermediate 
expressions, so they are not reported here. The hypothesis that (6.1) is a first integral 
uniquely leads to an equation that is equivalent to (4.7). 

The same method may be used to test whether our set of equations admits a constraint 
that would be a polynomial of degree 1 or 2 in the first derivatives. The only difference 
with respect to the procedure of looking for a first integral is that in verifying whether 
(6.2) is zero, eq. (6.1) is used, too. If a nontrivial solution of (6.2) with this additional 
simplification is found, then it means that the derivative of (6.1) by t is zero if (6.1) holds 
for any fixed t. Then, such (6.1) is a constraint preserved by the set S. However, even 
this attempt has not led to useful results. Constraints of degree 2, i.e. those with Qij ^ 0, 
lead to prohibitively complicated equations and could not be investigated. One constraint 
of the form (6.1) with Qij = was found, but it is equivalent to the square bracket in 
(4.11) being zero, and so implies zero matter density. Again, the details are not reported 
because they contain complicated equations, but no ingenious ideas. This result proves the 
usefulness of the method - a sensible constraint was revealed - but the solution with zero 
density is not interesting for cosmology, and thus not necessarily worth investigating. 

The zero-density constraint was found without using eq. (4.7). Eq. (4.7) would reduce 
the number of unknown functions by one, but the resulting set of equations is prohibitively 
complicated and no progress was achieved. 

7 Summary of results. 

It was shown that the Einstein equations for the metric (2.9) with / = g = are self- 
consistent and solvable. They reduce to the set S = {(4.3), (4.5), (4.6)} to determine h, 
K and K\\ = cosh(F), and (4.11) to determine ^33 (where C = 87rG/c 4 ). The matter 
density is found from (4.9). The first derivatives of the functions obeying the set S must 
obey (4.7). 

The Friedmann solution with k = — 1 is contained among the solutions of this set, as 
shown in eqs. (4.12) - (4.24). Unfortunately, no explicit example of a more general solution 
could be found. Attempts to follow ad hoc Ansatzes produced uninteresting results. The 
Ansatz K = K33 led to the deSitter solution in disguise, in which the t-lines had nonzero 




(6.3) 
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rotation. The Ansatz K\\ = K/C (C = const), which is consistent with the Friedmann 
limit, led to such complicated equations that it could not even be verified if they are 
not contradictory. The assumption of zero shear implies zero expansion, in virtue of the 
theorem (a = 0) (uj9 = 0) that holds for dust (see Ref. 12). 

The set S was shown to have a two-dimensional group of point symmetries, given by 
(5.14), and to admit no Lagrangian of the Hilbert type. It was also verified that no first 
integrals of the form (6.1) exist. 

The progress achieved in this paper was the reduction of the problem of existence of 
a rotating generalization of the k = — 1 Friedmann model to the technical problem of 
finding an explicit solution of the set S. The solvability of the set S may be taken for 
granted because the Friedmann model itself was shown to be one of its solutions. It is 
still unknown, though, whether a continuous family of solutions exists labeled by by the 
parameter uj (rotation) such that the limit u> — > taken in the explicit solution leads to 
the k = — 1 Friedmann model. 

A similar analysis as done here should be done for the other promising cases identified 
in Ref. 3. 

Acknowledgements. The algebraic manipulations for this paper were carried out 
using the computer algebra system Ortocartan 13 ' 14 that was for this purpose extended by 
several new programs 15 . The author is grateful to J. Kijowski, J. Jezierski, K. Rosquist, 
R. Jantzen, M. A. H. MacCallum and C. Uggla for useful comments and instructions on 
the Lagrangian methods. 

Appendix A 
The Einstein equations for the metric (2.9). 

As explained in sec. 2 (after eq. (2.9)), these are the projections of the Einstein 
tensor on the forms of the orthonormal tetrad (2.10), thus for example G03 below is equal 
to ei a ej^G a p with 1 = and J = 3, where G a p are the coordinate components of the 
Einstein tensor. 

Goo = 2K 11 - 3 hK 11 , t +K u - 3 K- 1 K 11 , t K, t +K 11 - 3 K 33 - 1 K 11 , t K 33 , t 
+^K 11 - 2 K~ 2 - ^K n - 2 K- 2 K 33 2 f, t 2 - ZK xx - 2 K- l hK, t 
-Ku- 2 K- l K 33 - l K, t K 33 , t -K xx - 2 K- x K, tt -^K u - 2 K 2 h, t 2 
-W u - 2 K 33 - l hK 33 , t -Kn^K^K^u +±K n - 2 K 33 2 hf, t g, t 
-±K n - 2 K 33 2 h 2 f, t 2 - X -K xx ~ 2 K 33 2 g, t 2 - 3A^rV - ^K u -%t 
+K n - 1 K- 1 K n , t K, t +K n - 1 K 33 - 1 K n , t K 33 , t - X -K~ 2 K 33 2 f, t 2 

+K- 1 K 33 - 1 K, t K 33 , t -3K- 2 (A.l) 
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Goi = -2K 1 r 2 hK 11 , t -K u - 2 K- 1 K u , t K, t -K 11 - 2 K 33 - 1 K 11 , t K 33 , t 
+^K 11 - 1 K- 2 + ^n-'r 2 ^ 2 // + K xl - l K- x hK, t 

+K ll - l K- 1 K, tt +K l r 1 K 33 ~ l hK 33 , t + K u - 1 K 33 - 1 K 33 , a +^K 11 ~ 1 h, t (A.2) 

—K^KK^K^t h, t ~K u - 2 Kh, tt ~K u - 2 K- 2 K, t 
+^K n - 2 K- 1 h + ^K^K^K^K^t +^K n - 2 K- 1 K 33 2 hf, t 2 
~K 11 - 2 K- 1 K 33 2 f, t g, t ~K u - 2 K, t h, t -K 11 - 1 K- 1 K 11 , t 

+2K~ 2 K, t -K-'K^K^t (A.3) 

G 03 = -^K n - 3 K 33 hK n , t f, t +^K n - 3 K 33 K n , t g, t +^K n - 2 K 33 hf, tt 
~K u - 2 K 33 hg, t +^K u - 2 K 33 h 2 f, t -±K n - 2 K 33 g, tt +^-K ll - 2 K 33 f \ t h, t 
+^K n - 2 hK 33 , t f, t -}-K n - 2 K- 2 K 33 f, t + l -K ll - 2 K- 1 K 33 hK, t f, t 

- X -K lx - 2 K- l K 33 K, t g, t -^K n - 2 K 33 , t g, t +^R- 2 K 33 f, t (AA) 

Gn = -K IX ' 2 K' 2 - jK u - 2 K- 2 K 33 2 f, t 2 + K xl - 2 K' x hK, t 
+K 11 - 2 K- 1 K 33 - 1 K, t K 33 , t +^K u - 2 K 2 h, t 2 + Ku~ 2 K^ 1 hK 33 , t 
~Kn- 2 K 33 2 hf, t g, t +^K u - 2 K 33 2 h 2 f, t 2 + ^K n - 2 K 33 2 g, t 2 
+K n - 2 h 2 - ±K u - 2 h,t -\K- 2 K 33 2 f, t 2 - K- l K 33 - l K, t K 33 , t 

-K~ l K, tt -Kaa^K^a +K~ 2 (A5) 
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G12 = -^Kn^KKn^Kt+^Kn^K-'Kn^+^Kn-'KK^K^th^ 
+±K n - 1 Kh, tt —Kn^K- 2 ^ +^if ! K' 1 K^ 1 K 33 , t 
—K^K^K^hfa 2 + X -K xx - x K- l K^f, t g, t ^-K xx - x K, t Kt (A6) 

Gl3 = ^ir 2 ^33^11,t/,t-^ir 2 ^33^11,t5,t-^ll _1 ^33/i/,« 

+^H- 1 J PC 3 35,« -^ii _1 ^33/,t /i,t -^ii" 1 ^,* /,t -Kn-'K-iRssU 

—Kn-'K^KsshK^ f, t +^K 11 ~ 1 K' 1 K 33 K , t g, t +^K n - 1 K 33 , t g, t (A.7) 

G 22 = -Kn^hK^t-Ku^Kas^K^tK^t+^Ku^K- 2 - ^K u ~ 2 K~ 2 K 33 2 f , t 2 

-jK n ~ 2 K 2 h, t 2 + 2K 11 - 2 K 33 - 1 hK 33 , t +K n - 2 K 33 - 1 K 33 , tt +±K n - 2 K 33 2 hf, t g, t 
-l Kll - 2 K 33 2 h 2 f, t 2 - jK n - 2 K 33 2 g, t 2 + K n ~ 2 h 2 + ^K n - 2 h, t 

-Kn-'K^K^t K 33 , t -Ku^Kum +\k~ 2 K 33 2 f , t 2 - K 33 - l K 33:tt +K~ 2 (A.8) 

G23 = ^K 11 - 3 K- 1 K 33 K 11 , t f, t +^K 11 - 2 KK 33 hf, t h, t ~K 11 - 2 KK 33 g, t h, t 

+^K n - 2 K' 2 K 33 K, t f, t -K xx - 2 K~ l K 33 hf, t 'K~' K 33 j \ u 

J-K lx - 2 K- l K 33 , t f, t +^K u - 1 K- 1 K 33 K u , t f, t -±R- 2 K 33 K, t f, t 

+^K' 1 K 33 f, tt +^K~ 1 K 33 , t f, t (A.9) 

G 33 = -Kn^hK^t -K^K^K^tK* --K^K' 2 + \k u ~ 2 K~ 2 K 33 2 f \ t 2 
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+2K 11 - 2 K- 1 hK, t +K 11 - 2 K- 1 K, tt —^K\\~ 2 K 2 h, 2 + ^K u - 2 K 33 2 hf, t g, t 

—Ku^K^h 2 /,^ - -K xl - 2 K Z3 2 g, t 2 + K u ~ 2 h 2 + ^K u ' 2 h, t -^f 1 K' 1 K u , t K, t 

-Ku^K^u -^K- 2 K 33 2 f, t 2 - K~ X K M +K~ 2 (A10) 

Since the source in the Einstein equations is dust with a cosmological constant, and since 
the zero-th tetrad vector is just the velocity vector, the above components should obey the 
following equations: 

Goo = (8vrG/c 4 )p - A, 

Gii = G22 = G33 = A, 

nondiagonal G/j = 0, 
where p is the dust energy-density and A is the cosmological constant. 

Appendix B 

Consequences of h = in the Einstein equations. 

With h = 0, eq. (2.12) becomes: 

K 33 f,t (K n 2 ~ l)/(K n K 2 ) = (B.l) 

We can immediately discard the solution K 33 = because then det{g a p) = 0. When 
if ii 2 = 1, the limit oj — > of the resulting metric will necessarily have either nonzero shear 
or zero expansion (see sec. 3 for a calculation of this limit), and so no generalization of 
the Friedmann models can be expected here. Hence, the only consequence of (B.l) that is 
worth pursuing is: 

f,t = 0. (B.2) 

Then, with h = 0, eq. (A. 6) implies: 

KuK 33 /K 3 = const. (B.3) 

However, in the limit oj — > this again implies either nonzero shear or zero expansion, i.e. 
no Friedmann limit. 

Appendix C 
The result for F, tt from Gqi = 0. 
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When (4.2), (4.2) and (4.3) are substituted in (A. 2), the following formula results for 
F,tt- 



F, tt = ~Khcosh- 1 (F)smh- 1 (F)K, t h, t +3hcosh' 2 (F)F, t - 7 -hsmYr 2 {F)F, t -3hF, t 

+ ^K- 2 cosh- 1 (F) sinh(F) + 2K- 2 cosh- 1 (F) smh(F)K, t 2 
+^K- 1 hcosh- 1 (F)smh~ 1 (F)K, t +3K- l hcosh~ 1 (F) smh(F)K, t 
+3K~ 1 cosh~ 2 (F)K, t F, t +3K- 1 cosh~\F) smh(F)K, tt -lK~ x K,t F, t 
+^K 2 hsmh- 2 (F)F, t h, t -3K 2 h 2 cosh- 1 (F)smh' 3 (F)h, t 
--K 2 h 2 cosh- 1 (F)smh- 1 (F)h,t--K 2 cosh-\F)sm-h- 1 (F)h, t 2 



+/i 2 cosh" 1 (F)sinh" 3 (F) + -h 2 coBh- 1 {F)gmir 1 (F) + 2cosh- 1 (F)sinh- 1 (F)F,t 2 



1 3 

+-cosh -1 (F)sinh -1 (F)/i,t +2cosh -1 (.F) sinh(F)F, t 2 + -cosh _1 (F) wah(F)h, t . 

It will be modified later because it contains K,tt on the right-hand side, while the final 
equations that will be dealt with should have no second derivatives on the right-hand 
sides. 

Appendix D 

Nonexistence of a Hilbert-type Lagrangian for the set {(4.3), (4.5), (4.6)}. 

Eqs. (4.3), (4.5) and (4.6) can be written in the form: 

d 2 z l TTr , dz J dz k ■ dz j TTj 

where i = 0, 1, 2; z° = h, z 1 = K, z 2 = F and W % jk, V l j and U l are functions of (h, K, F) 
(but not of t). Note that the set (D.l) is covariant with respect to arbitrary transformations 
z % — > z H = f l ({z^}): the first derivatives ^jj- transform then like a contravariant vector, and 
so do the terms U l , the coefficients V l j transform like a mixed tensor, and the coefficients 
(— W l jk) transform like components of an affine connection. (The nontensorial terms in 
the transformed (—W l jk) arise from ^p-)- The most natural ansatz for a lagrangian for 
(D.l) is: 

„ dz 1 dz 3 r dz l 

L = Qu- — + Li— + <5>, (D.2) 

^ 1 dt dt dt v ' 
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where Qij, Lj and <I> are functions of (h,K,F). Such a lagrangian would result from the 
Hilbert lagrangian by taking out a complete divergence and integrating the result with 
respect to the spatial variables. The Euler-Lagrange equations implied by (D.2) are: 

d 2 z s 1 dz k dz l 1 dz k 1 /r»o\ 

^ is "dP" = ~ 2^'^"dT"dT + 2^ ~~ ' fc '~d7 + 2 ^ ^ 

If these are to be equivalent to (D.l), then the following must hold: 

QisW s M = ~(Qki,i + Qh,k ~ Qkl,i), (DA) 
QisV S k = -^{Lk,i - L it k): (D-5) 

QisU s = \$,i. (D.6) 

Eqs. (D.4) imply that (— W l jk) must be Christoffel symbols constructed from the metric 
Qij, eqs. (D.5) imply that \Li must be a vector potential for the tensor field Qi s V s k, and 
eqs. (D.6) imply that <J>/2 must be a scalar potential for the vector field Qi s U s . All of 
these are strong conditions and they may be impossible to fulfil in many cases. 

Indeed, for our eqs. {(4.3), (4.5), (4.6)}, the solution of (D.4) turns out to be Qij = 0, 
i.e. the Lagrangian (D.2) does not exist. This is an outline of the proof. 

After eqs. (D.4) are written out in the form 

Qij,k = -W s ik Q S] - W s jk Q ls , (D.l) 

with W s ki read off from {(4.3), (4.5), (4.6)}, the following two equations follow, among 
other results: 

Q lltF + ^(K cosh(F)/ smh(F))Q n ,K = -(2cosh 2 (F) - l)Q n /(cosh(F) sinh(F)), (D.8) 



Q22,K + (2cosh 2 (F) - l)sinh(F)Q 22 ,F/(2^cosh 3 (F)) = (3cosh 2 (F) - l)Q 22 /(Kcosh 2 (F)). 

(D.9) 

The solutions of these are: 

( K A \ 1 
Qii=qii[h,——\ . (D.10) 

\ smn(f ) J cosh(i < ) sinh(i< ) 



/ KJ2cosh 2 (F) - l\ 
Q 22 = Ksmtf(F) q22 k V s . nh2(F) J , (D.ll) 

where are arbitrary functions of their two arguments. The equation Qu t K = ... is then 
solved with the result: 

Q12 = -K 5 q lliW / S mh. 3 (F), (D.12) 

where w = K A / sinh(F) is the second argument of qn, and the equation Q 22) k = ... 
implies: 
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^22,« = qii, w K 3 /(cosh(F) S mh 2 (F)), (D.13) 

where v is the second argument of g22- The left-hand side of (D.13) is an invariant of 
the operator (2Kcosh 3 (F)/ sinh(i ? ))^ + (2cosh 2 (i ? ) — 1)^, while qii tW is an invariant of 
the operator \K-^ + (sinh(F)/ cosh(i ? ))^. Application of these two operators to (D.13) 
leads to (722,1; = Q.u,w = 0, which implies Q12 = 0. With this, the remaining equations 
(D.7) quickly lead to Qij = 0, which means that the Lagrangian (D.2) does not exist in 
this case. 

Since the Euler-Lagrange equations (D.4) are covariant with respect to arbitrary trans- 
formations of the Lagrangian variables (in our case h — > h'(h, K, F), etc.), and equa- 
tions of the form (D.l) are covariant, too, the conclusion that a Lagrangian of the form 
(D.2) exists (or does not exist) is coordinate-independent, i.e. having shown that eqs. 
{(4.3), (4.5), (4.6)} do not follow from a Lagrangian (D.2) in our variables {h,K,F}, we 
know that no such Lagrangian will exist in any other variables. 

Appendix E 
The general solution of eqs. (5.8) — (5.11). 

Eqs. (5.8) have the form: 

t-M = 0, (E.l) 

where ; is the covariant derivative in which (— W\j) play the role of the connection co- 
efficients. (They appear in this role for a second time already, see Appendix D.) The 
integrability conditions of (E.l) are: 

R'ijk^s = 0, (E.2) 

where R s ijk = —R s ikj is the curvature tensor corresponding to the connection(— W 1 ^). 
Eqs. (E.2) are 9 equations (labelled by the sets of indices k) = (0, 0, 1); (0, 0, 2); (0, 1,2); 
etc) and they could have nontrivial solutions only if every subset of 3 equations chosen 
from among them had a zero determinant. Actually, of the 84 determinants only two 
vanish, and some of them will not vanish even if the functions h(t), Kit) and F(t) are 
functionally dependent. Here is one determinant that will never vanish, it corresponds to 
{(i,j,k)} = {(1,0,1); (1,1, 2); (2, 1,2)}: 

det(£.2) = K-\- 1 ^- + ^cosh- 6 (F) + ^cosh- 4 (F) + ^cosh- 2 (F) - ^sinh^F))]. 
Hence, the unique solution of (5.8) is: 

With £,j = 0, eqs. (5.9) simplify somewhat, and the equation corresponding to (i, k, I) = 
(0, 1,2) becomes t]°,kf= — 2??°,f /K, which has the solution: 

rf = F°(t,h,F)/K 2 + G°(t,h,K), (EA) 

F° and G° being unknown functions. Then, eq. (5.9) with (i,k,l) = (0,1,1) allows us to 
separate the variables K and F, and its solution, substituted into (E.4), gives the result: 
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rf = M°(t, h)smh\F)/K 2 + J°(t, h)/K 4 + L°(t, /i), (E.5) 

where M , J° and L° are new unknown functions. With this, eq. (5.9) corresponding to 
(i, k, I) = (0, 2, 2) implies J° = 0, and the one with (i, k, I) = (0, 0, 2) solves as follows: 



„> . ^Mo g(8 inh(F))- Am% ^HF) + ^ f I2; { ^ ^ *■ *)■ <*«) 

where /i, K) is a new unknown function, and ?7 2 is still completely unknown. Then, 
for (i,k,l) = (0,0, 1), the equation (5.9) has the solution: 

F 1 = -^M°KhlogK + G 1 (t,h)K, (E.7) 

where G l {t, h) is a new unknown function. 

When (E.6) and (E.7) are substituted into the (0, 0, 0) component of (5.9), an algebraic 
equation for rf results, whose solution is: 



rf = i (-^^cosh(F)sinh 5 (F) _ ^^^cosh(F)sinh 3 (F) 

' 3cosh 2 (F) + l \3 K A h y ' y ' 3 K 2 h V ' K ' 

35 M° 5 L° 

H ^- cosh(F) sinh 3 (F) H — cosh(F) sinh(F) 

6 K z h 2 h 

+3M°/i cosh(F) sinh(F) [log(sinh(F)) - log K] + 5 cosh(F) sinh(F)G 1 
cosh(F) sinh 3 (F) + h°, h cosh(F) sinh(F) 



" 5 S cosh ( F ) sinh 3 (F)[log(sinh(F)) - log if] - y^f cosh(F) sinh 3 (F) . (E.8) 

Both sides of the (1, 1, 1) component of (5.9) become then polynomials in (log if) and 
1/K, whose corresponding coefficients have to be respectively equal. The coefficients of 
K^ 1 logiC imply then M° = 0, and with this, only two other terms remain whose solutions 
are: 

L° = 3C(i)/i~ 5/3 - B(t)h, G 1 = C(t)h~ 8 / 3 + B(t). (E.9) 

where C(t) and B(t) are unknown functions. But this results in rf = 0, rp = L°, rf = G 1 K. 
Then, the (1, 0, 0) component of (5.9) implies C(t) = 0, and with this all the remaining 
equations in (5.9) are fulfilled. Thus the final solution of (5.9) is: 

v ° = -B(t)h, rf = B(t)K, rf = 0. (£7.10) 
With £ = £(t) from (E.3) and rf as above, any equation of the set (5.10) implies: 

B = const, £ = Bt + A, A = const, (E.ll) 
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and this satisfies all the remaining equations (5.10) and (5.11). Hence, the general solution 
of (5.8) - (5.11) is (5.12). 

This result was derived under the tacit assumption that the functions h, K and F are 
functionally independent. In the course of solving the equations (4.3), (4.5) and (4.6), 
relations between these functions may appear. It happens sometimes that such relations 
are revealed by the symmetry equations as cases in which the symmetry group is larger 
than in the generic case (see e.g. Ref. 16 where special cases of larger symmetry of a 
single equation were revealed by the symmetry equations). This possibility has not been 
investigated for the equations (5.9) - (5.11). However, for the equation (5.8) the solution 
is always (E.3), even if the functions h, K and F are not independent, as shown in the 
paragraph containing (E.3). 
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